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Analytical Result on Electronic States around a Vortex Core in a 
Noncentrosymmetric Superconductor 
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We analytically study electronic bound states around a single vortex core in a noncentrosym- 
metric superconductor such as CePtsSi without mirror symmetry about the ah plane. Consid- 
ering a mixed spin-singlet-triplet Cooper pairing model, we obtain a formula about the local 
density of states (LDOS) around a vortex core in any direction of the magnetic field. The LDOS 
under a magnetic field perpendicular to the c axis is quite different from that of typical s-wave 
or d-wave superconductors. From the ellipticity of the spatial pattern of the LDOS around a 
vortex core, one can experimentally estimate the pairing symmetry of CePtsSi, such as the 
position of the gap nodes and the ratio of the singlet component to the triplet component in 
the order parameter. 

KEYWORDS: CePtaSi, Unconventional superconductivity, Vortex core, Local density of states, Broken 
inversion symmetry, Quasiclassical theory 
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In a time reversal invariant system, the lack of in- 
version symmetry is connected with the presence of 
an antisymmetric spin-orbit coupling. -'^ Unconventional 
superconductivity in strongly correlated systems exists 
in heavy Fermion superconductors. The discovery of 
the noncentrosymmetric heavy Fermion superconduc- 
tor CePtaSi^ has attracted considerable attention. Many 
theoretical and experimental studies on unconventional 
superconductivity without inversion symmetry have been 
reported for a few years. The vortex structure of the 
mixed state in this system was recently studied on the 
basis of the Ginzburg-Landau theory and the London 
theory.^'"' The vortex core structure was also studied 
by means of the quasiclassical theory of superconduc- 
tivity, which enables us to calculate the physical quan- 
tities more microscopically.^ However, no studies have 
been attempted to investigate the possible anisotropic 
spatial structure of the quasiparticle distribution around 
a vortex core in such a system. In this paper, we ana- 
lytically investigate the local density of states (LDOS) 
around a single vortex core in a noncentrosymmetric su- 
perconductor and devise a formula for LDOS that is ap- 
plicable in any arbitrary direction of the magnetic field, 
thus revealing an anisotropic quasiparticle structure. An 
anisotropic pairing symmetry is reflected by the real- 
space LDOS pattern because of the presence of quasi- 
particles around a vortex core.^^'^^ Therefore, for un- 
derstanding the Cooper pairing without inversion sym- 
metry, it is important to investigate the LDOS pattern. 
The LDOS can be probed experimentally by STM.i3^i4 
If the LDOS pattern presented in this paper is consistent 
with that observed by STM, we can obtain information 
on the position of the gap nodes and the ratio of the sin- 
glet component to the triplet component in this material. 
In other words; this can serve as one of the methods to 
experimentally investigate the pair potential directly in 



CePtaSi. Thus, we can also show that our theoretical 
formulation is capable of wide application. 

We consider a mixed spin-singlet-triplet model to 
study the noncentrosymmetric superconductor. Consid- 
ering the results of numerical calculations by Hayashi 
et al.,^ we assume that the spatial variations of the s- 
wave pairing component of the pair potential are the 
same as thos of the p-wave pairing component. A = 



hy d'x 



kx d'y 



[^ao + dk-&] i&y ^ A(r) 
where the s-wave pairing component 4*, the d vector 
dk = A(— fcj,, kx, 0), the Pauli matrices in the spin space 
o' = {d'x, cTy,az), the unit matrix (Tq, and the unit vector 
k in the momentum space will be discussed later. Here, 
the ratio of the singlet to the triplet component is defined 
in the bulk region as 4* = 'f/A. This mixed s-|-p-wave 
model is proposed for CePtsSi.^ 

In a system without inversion symmetry, there is a 
Rashba-type spin-orbit coupling with the form"*' 



Hi 



k,T],r]' 



(1) 



where flfe = y fcx, 0). Here, c^(cfc^) is the cre- 

ation (annihilation) operator for the quasiparticle state 
with momentum k and spin rj and a (>0) denotes the 
strength of the spin-orbit coupling. We use units in which 
h ^ ks = 1. The gTj, vector, which is an antisymmetric 
vector = —Qk), is parallel to the d vector.* 

We calculate the LDOS around a vortex core on the ba- 
sis of the quasiclassical theory of superconductivity. ^^"^^ 
We consider the quasiclassical Green function g that has 
matrix elements in the Nambu (particle-hole) space as 



g{r,k,iujn) = 



9 I 
-/ 9 



(2) 



where ujn is the Matsubara frequency. Throughout the 
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paper, "/lai" a denotes a 2 x 2 matrix in the spin space, 
and check" a denotes a 4 x 4 matrix composed of the 
2x2 Nambu space and the 2x2 spin space. 

The Eilenberger equation, which includes the spin- 
orbit coupHng term, is given as^' i^^^o 



ivF{k) ■ Vg + [iuinfs - ag^ ■ S - A, g] = 0, (3) 



where 



9k = 



_ ( 9k^o 



g-k^o 



9k^o 
-9k^o 



9k = V 2(~^2/'^^'0)' ^3 



(T 





A 



ao 
-ao 

A 
-At 



, (4) 



(5) 



(6) 



ay&ay, and 



Here, VF{k) is the Fermi velocity, o-*"^ : 
the commutator [d, 6] = ah — ha. g should be subject to 
the normalization condition^^' g^ = 1, where 1 is a 
4x4 unit matrix. We neglect the impurity effect and 
the vector potential because CePtsSi is a clean extreme- 
type-II superconductor.^ 

The Eilenberger equation can be simplified by intro- 
ducing a parametrization for the propagators that satisfy 
the normalization condition. Propagators are defined as 
P± = i (i =p g), which were originally introduced in the 
studies of vortex dynamics. Using these propagators, 
we obtain the matrix Riccati equations as follows: 



Vf 



■ Va+ -|- 2u)na+ + a+A'^d+ — A 



+ i{a+a{g^ ■ af' + ag,, ■ &a+) 
Vf ■ VS_ - 2u;„&_ - S_A6_ + A^ 



0, (7) 



0, (8) 



where 

9 = 

N 



= -N 



(i - a+L 
-2ib- 



(1 + a+h. 




2ia+ 
-(i - 6_a+) 



(l + 6_a+)-i 



, (9) 



(10) 



We consider a single vortex along the Z axis that tilts 
from the crystal c axis by an angle (j). Now, we con- 
sider the X axis on the a-b plane. We assume the spher- 
ical Fermi surface and consider the momentum vector 
k = A;(cos0sinx, sin6'sinx, cosx) in this coordinate sys- 
tem fixed to the magnetic field. To obtain the quasi- 
classical Green functions, we consider the trajectories 
of the quasiparticle on the X-Y plane. ^'^ Because of a 
translational symmetry along the vortex, the trajectory 
of a quasiparticle with a momentum |A;f| and kz ^ 
contributes to the quasiclassical Green functions in the 
same manner as that of a quasiparticle with a momentum 
jfepl sinx projected on the kx-ky plane. In other words, 
the quasiclassical Green functions in a three-dimensional 
system can be converted into a set of quasiclassical Green 
functions in a two-dimensional system having a momen- 
tum with different amplitude. Therefore, we determine 



the coordinates 

X 
Y 



cos u 
sin 6 



— sin 9 
cos 6 



(11) 



to solve eqs. (7) and (8) along the direction of momentum 
in a two-dimensional system; the momuntum is in the x 
direction. Here, y is the axis referred to as an impact 
p arameter , x is the axis perpendicular to y, and r = 
■\/x^ +y'^ = \/X'^ -h y2. ^i^g origin of these axes is at the 
vortex center. 

Using the perturbative method developed by Kramer 
et aL,^^~^* we can obtain the quasiclassical Green func- 
tion around the vortex core in a low energy region 
(|a;„| <C |A(X)|). Here, Aqo is a pair potential in the bulk 
region. By expanding the matrix Riccati equations in the 
first order of a, y, | A(r)|, and |a;„|, we obtain the approx- 
imate solution as^^ 



Q^^-1{^k^k.+^)F{x) 



2F. 



+ 



->S>\F{x) 



2Fc 



(12) 



where 

Fc,D 



1 /"^ 
- / dx' 

|A(r' 



)l y 



X e 



smx 

-2j^fc+fc_^*|_F(a;' 



± 



C,D 



(13) 



(14) 



Here, k± = ka,± iky and F{x) = J^^' dx'A{x')/vF. The 
first and second terms in cq. (12) arc the Green functions 
referred to as gi^n on Fermi surfaces I and II; this is 
because the spin orbit coupling splits the Fermi surface 
into two surfaces by lifting the spin degeneracy.^' '"^ Fermi 
I has a nodeless pair potential, and Fermi II has a pair 
potential with the line node (e.g., see Fig. 1 in ref. 10 ). 

Near the vortex core (|r| <C ^), we use the following 
approximations. 



Fc,D 
where 



y\^k+k-T'^\ 
sinx 



+ 



^,(15) 



F{x) 



-f 

Vf J-< 



- /' \x'A(x') ~ 0, 
Vf Jo 



Jo 



dx' 



aAoo 

Acx,. 



(16) 
(17) 
(18) 



Here, we use the dimensionless variables as d)„ = Wn/ A^o 
and y = y/^o with = wf/Aqo. Therefore, the LDOS 
around the vortex c;orc is 



r>(f,e) = -(Re{trff^'=(jw„^? + 0+)}>. 



(19) 



z>i(0) 



sinxdxj^ — (x/fc+fc- +*) 
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xs[h+ie,x)] 

+ i>ii(0) / sinxdx 

Jo Jo 



\/k+k- 



x6[h^i0,x)], (20) 

where x) = e — y{\/k+k^±^)^ /sinx- Here, t'i,ii(0) 
denotes the density of states on Fermi surfaces I and II. 
This integral reduces to an integral on the path satisfy- 
ing h±{9,x) = 0. Therefore, at the points where v{f,e) 
diverges, the following relations hold: 

dh± dh± 



h± = 







dO dx 
and the saddle point condition 



<0. 



(21) 



(22) 



a6»2 dx^ 

These relations are derived from the following discussion. 
Near the saddle point, z>(f, e) is rewritten as 



JV3{0,0) 

with ab < 0, and diverges logarithmically. 
From eq. (21), it follows that 



(23) 



x± = 



y± = 



= 



2esmx-§g^/W+H 



fc2 + fc2 ± * 



esmx 



^Jkl + kl±^ 



ecosx 
sinx 



-dx 



{kx + ky) 



i.2 



i.2 



1.2 



(24) 



(25) 



(26) 



We can consider eq. (26) as the equation to determine 
X- The solution of eq. (21) is regarded as an envelop- 
ing curve of the trajectories of the quasiparticle when x 
is fixed. 2^ Therefore, from eq. (11), we can obtain the 
points where the LDOS diverges. 

In momentum space, the relation between the axes 
fixed to the crystal axes {kx,ky,kz) and those fixed to 
the magnetic field {kxjky, kz) is written as 



kx = kx, ky = ky cos ^ — kz sin (f 
kz = ky sin (j) + kz cos ( 



(27) 
(28) 



Here, we consider the rotation about the a axis; this ro- 
tation reveals that the X axis is equal to the a axis. 
Therefore, using the spherical coodinates fixed to the 
magnetic field {kx = A;f cos ^ sinx, ky = fcp sin ^ sinx, 
kz = kp cosx), we obtain 



kx 
ki, 



cos ti sm X, 
sin 9 sin x cos ( 



cos X sm ( 



(29) 
(30) 



First, we consider the system in a magnetic field par- 
allel to the c axis (^ = 0). Equation (26) leads to x = f • 



Therefore, we obtain the parametric equations 



esinO 
'(1±*): 



and Y± 



ecos0 

(1±^')2 



(31) 



from eqs. (22), (24), and (25). This result is consistent 
with the numerical calculation by Hayashi et al.^ because 
this spatial LDOS pattern is in the form of two concentric 
circles. This result shows that the quasiparticles rotate 
around a vortex core, and the LDOS exhibits the two- 
gap property. The ratio of the singlet component to the 
triplet component, determines the ratio of the radius 
between the two concentric circles. 

Second, we consider the system in a magnetic field 
perpendicular to the c axis (</> = ^). In this case, X = f 
follows from eq. (26). Therefore, we obtain the following 
parametric equations: 



— esmf 



(|cos6'| ±*)3 



|3|cos6'|±^'|, 



e cos 6 



^ 2 ^^^^ 6'sign[cos 6*] 

(|cos6'| ±*)3 ^ {\cose\ ±*)2' 



(32) 
(33) 



where ^E* < |cos0| ( in the case of /i_) using eqs. (22), 
(24), and (25) with regard to Fermi I {h+) and II (hJ). 
In Fig. 1, we show these LDOS patterns aroimd the vor- 
tex core for a fixed energy e. The pattern about Fermi 
I in Fig. 1(a) can be regarded as a result of the shape 
of the order parameter, such as two circles with a shift 
of the centers (shown in Fig. 1(a) in ref. 10). If the re- 
lation ^ < I cos ^ I between ^ and 6 does not exist in 
the case of the pattern about Fermi II is similar 
to that in a d-wave superconductor^^' ^2 because of the 
presence of four nodes in the pair potential on the cir- 
cular line X = ■! on the Fermi surface. However, the 
actual pattern about Fermi II (Fig. 1(b)) is similar to 
a part of that in a d-wave superconductor because of 
the presence of this relation between and 9 originat- 
ing from the three-dimensional anisotropic pair poten- 
tial. These results show that the quasiparticles on Fermi 
I are bound aroimd a vortex, while those on Fermi II 
move in the direction associated with one node to the 
direction associated with the other node. We also cal- 




Fig. 1. Distribution of the points where the local density of stctcs 
diverges around the vortex perpendicular to the c axis, (a) Fermi 
I {h+) and (b) Fermi II (h-). Here, * = 0.5 and e/Aoo = 0.01. 
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culate the distribution of the LDOS from eq. (20) by 
numerical integration (Fig. 2). The LDOS in Fig. 1 is 
consistent with that in Fig. 2. In the Umit ^' — > 0, the 
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Fig. 2. Distribution of the LDOS around the vortex perpendicu- 
lar to the c axis, (a) Fermi I (/»+) and (b) Fermi II (h-). Here, 
'I' = 0.5 and e/Aoo = 0.01. Notice the diflference in spatial extent 
between these LDOS patterns. 



LDOS patterns in a magnetic field perpendicular to the 
c axis become rather isotropic because the line nodes in a 
gap disappear. Therefore, the strongly anisotropic LDOS 
patterns, as shown in Fig. 2(b), suggest that the triplet 
channel and the singlet channel are mixed. 

We can estimate the ratio of the singlet component to 
the triplet component, ^P, from the spatial LDOS pattern 
shown in Fig. 2(b). From eqs. (32) and (33), we calculate 
the ratio of the intercept Xq on the X axis to the inter- 
cept Yq on the Y axis in the LDOS pattern for Fermi II. 
The ratio between these intercepts r is written as 



r = 



(-5^ + V2A + l^'f 



108(1 - ^)2(-§ + + 24)^1 - (^*±V|!±24 

(34) 

where r = Yq/Xq (Fig. 3); this ratio does not depend 
on the energy e. In other words, we can obtain the ratio 
from the cUipticity r of the shape around a vortex. 
We can also obtain the position of the gap-node ^node 
with the relation |cos^node| = ^- Now, CePtsSi has a 
tetragonal crystal structure. Wo have also investigated 
how the present result is influenced by an anisotropic 
mass tensor; thus, we found that the LDOS pattern is 
not significantly affected by the ellipsoidal Fermi surface. 

In conclusion, we investigated the local density of 
states around a vortex core in a nonccntrosymmctric su- 
perconductor. We derived an analytical formula of the 
LDOS in any direction of the magnetic field. In a mag- 
netic field parallel to the c axis, we found that the LDOS 
pattern is in the form of two concentric circles; this is 
consistent with the numerical calculation by Hayashi et 
al.^ In a magnetic field perpendicular to the c axis, we 
found that the LDOS pattern about Fermi I is as shown 
in Fig. 1(a) and that about Fermi II extends far from 
the vortex center as shown in Figs. 1(b) and 2(b); this 
pattern is similar to that of a d-wave superconductor, 
but it is different from a four-fold d-wave superconduc- 
tor. The anisotropic LDOS patterns indicate the mixed 
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Fig. 3. Relation between the ellipticity r of the pattern around a 
vortex in the LDOS pattern about Fermi II and the ratio of the 
singlet component to the triplet component, ^. 



singlet-triplet channels. We can obtain the ratio of the 
singlet component to the triplet component in CePtaSi 
from the ellipticity of the shape around a vortex as shown 
in Fig. 1 (b) . While we found that the effect of the defor- 
mation from the spherical to the ellipsoidal Fermi surface 
is not quite significant, it may be important in the fu- 
ture to investigate a more realistic band structure. Our 
analytical formulation presented in this paper is advan- 
tageous because it can be easily generalized to a system 
with anisotropic Fermi surfaces. 
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